Rotonlike instability and pattern formation in spinor Bose-Einstein condensates 
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We show that metastable phases of an antiferromagnetic spin-1 condensate in a simple model with 
pure contact interactions can exhibit a rotonlike minimum in the excitation spectrum. The intro- 
duction of magnetic held gives rise to the instability of roton modes, which can lead to spontaneous 
emergence of regular periodic, polygonal, polyhedral or crystalline patterns, as shown in numeri- 
cal simulations within the truncated Wigner approximation. An explanation of the occurrence of 
rotonlike instability is given based on the energy and spin conservation laws. 
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Bose-Einstein condensates with spin degrees of free- 
dom [H attracted in recent years great interest due 
to the unique possibility of exploring fundamental con- 
cepts of quantum mechanics in a remarkably controllable 
and tunable environment. The ability to generate spin 
squeezing and entanglement 0] makes spinor Bose gases 
promisiiig candidates for applications as quantum sim- 
ulators 0, in quantum information and for precise 
measurements [5|- Moreover, spinor condensates were 
successfully used to recreate many of the phenomena of 
condensed matter physics in experiments displaying an 
unprecedented level of control over the quantum system. 
In particular, spin domains spin mixing Q, and spin 
vortices |8i] were predicted and observed. 

The fundamental concept of a roton excitation, first 
introduced by Landau in the context of superfluid he- 
lium, is crucial for understanding of its physical prop- 
erties It is characterized by a minimum in the 
spectrum of excitations occurring at a finite wavelength, 
E{k) ss A -I- h?{k — fco)^/2/Lt. If the roton gap A can be 
decreased by changing the system parameters, the soften- 
ing of the roton mode can eventually lead to an instabil- 
ity. This instability scenario is encountered also in many 
other branches of quantum physics, including stro ngly 
correlated Fermions [13], quantum Hall systems 11 1, 
and Bose-Einstein condensates with long range interac- 
tions 12h14|. The roton instability is characterized by 



unstable modes with wavevector lengths close to the ro- 
ton minimum fco- It was suggested that it can lead to 
the emergence of the peculiar supersolid state [15|, and 
several other physical phenomena 16 1 . 

Here, we show that rotonlike instability can occur in 
spinor Bose-Einstein condensates in a simple model with 
pure contact interactions. It can take place in appropri- 
ately prepared metastable states of an antiferromagnetic 
spin-1 condensate 17] under the influence of magnetic 
field. Moreover, we show that, depending on the geome- 
try and the trapping potential, it can lead to spontaneous 
emergence of variety of transient ordered patterns, in- 
cluding polygonal, polyhedral and crystalline structures. 
We show that these results can be verified in experi- 



ments with ^^Na condensate, which is characterized by 
very weak dipolar interactions [l^]- We provide an ex- 
planation for the occurrence of rotonlike instability based 
on energy and spin conservation laws, and demonstrate 
how the pattern characteristic length is determined by 
the transfer of the spin energy to the kinetic energy. 

We consider a dilute spin-1 BEC in a homogeneous 
magnetic field pointing along the z axis. We start with 
the Hamiltonian H = Hq + H^, where the symmetric 
(spin independent) part is 



i=-,o,+ 



(1) 

where the subscripts j = — ,0,+ denote sublevels with 
magnetic quantum numbers along the z axis m/ = 



— 1,0;+1, m is the atomic mass, n 
is the total atom density and V{r) - 



\z^ is the external potential. The asymmetric part 



can be written as 



dr 
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where Ej are the Zeeman energy levels and the spin 
density is F = {-^^ Fx'il>,'il)^ Fyip,'4)^ Fz'4>) where F^^y^^ are 
the spin-1 matrices and ip = ('!/'_|_, V'o, V'-)- The spin- 
independent and spin-dependent interaction coefficients 
are given by cq = 47rfi.^(2a2 -|-ao)/3m and C2 = 47r/i^(a2 — 
ao)/3m, where as is the s-wave scattering length for col- 
liding atoms with total spin S. 

The linear part of the Zeeman effect induces a ho- 
mogeneous rotation of the spin vector around the di- 
rection of the magnetic field. Since the Hamiltonian is 
invariant with respect to rotations around the z axis, 
we can remove this trivial effect by introducing a frame 
rotating with the Larmor frequency. We thus consider 
only the effects of the quadratic Zeeman shift For 
sufficiently weak magnetic field we can approximate it 
by a positive energy shift of the m/ — ±1 sublevels 
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FIG. 1: (a) Energy spectrum of stationary states of a uni- 
form spin-1 antiferromagnetic condensate with magnetization 
M = as a function of the magnetic field strength. The scal- 
ing field Bo depends on the density, (b) The phase-matched 
(PM) phase is an excited state, however in the case _B = it 
corresponds to maximal transverse polarization and its stabil- 
ity is ensured by spin conservation. For B > the transverse 
spin is no longer conserved and instability occurs, (c) Exci- 
tation spectrum of the PM phase in the case B — 0. The 
dash-dotted, dashed, and solid lines correspond to density 
wave, spin wave and spin quadrupole Bogoliubov modes, (d) 
The same for B = l.lGBo. The minimum in the quadrupole 
mode gives rise to an instability characterized by imaginary 
frequencies. The spin healing length is = 2nh/y'2mc2n. 



SE = {E+ + E_ - 2Eq)/2 « a2£:HFs/16, where £:hfs 
is the hyperfine energy splitting at zero magnetic field, 
a — {gi + gj)nBB/EYiFS, A*s is the Bohr magneton, gj 
and gj are the gyromagnetic ratios of electron and nu- 
cleus, and B is the magnetic field strength 17 [. 

In the mean field approximation, the Hamiltonian gives 
rise to the Gross-Pitaevskii equations 



= ['C + C2{n± + no - n=p)] ■ip± + C2ipo'>P^ , (3) 

= ['C - (5£^ + C2(n+ + ?!_)] V'o + 2c2^p+'lp-i>o ' 

where C = —h^V^ /2m + cqu + V{r). The total number 
of atoms N — J ndr and the total magnetization AI = 
/ Fzdr — J (n+ — n_) dr are constants of motion. 

First, we describe possible phases of a condensate in 
a homogeneous magnetic field in the case of a vanish- 
ing potential, V(r) = [13] ■ Homogeneous stationary 

■g-iAij t-i-i0j ^ where 



solutions have the form ij}j{Y,t) - 
6j are the phase shifts and /i4. -I- /i_ — 2/io- These so- 
lutions are stationary in the sense that the number of 
atoms in each magnetic sublevel rij is constant in time, 
but the relative phases may change as a result of an 



additional spin precession around z. In Fig. [TJa) we 
show the energy per atom in various phases of an an- 
tiferromagnetic condensate (c2 > 0) as a function of 
the magnetic field strength for M — 0. Because the 
symmetric part of the Hamiltonian ([1]) is a constant, 
the relevant part of the energy e(r) is defined through 
Ha = Jdrne{r). Note that due to scaling properties 
of the Hamiltonian, the use of renormalized variables 
e/{c2n) and B/Bq, where Bq depends on the density 
through the condition SE\b=Bo = C2n, allows to show 
the possible phases of spin-1 condensate in one universal 
graph (there is no fixed parameter except M). The ne- 
matic (po) state is described by -0 = (0, y/n,0), the two- 
component (2C) state hy ip = {y^n/2, 0, A/n/2), while in 
the phase-matched (PM) state all three magnetic com- 
ponents are populated with the relative phase between 
them equal to 9 = 9+ + 6- — 29^ — and relative popu- 
lations dependent on the magnetization. In the limit of 
zero magnetic field, nematic and 2C states become polar 
states, while PM states become ferromagnetic states |l|. 
For more details about the possible phases, see [l7| . 

The PM phase is the highest excited state at low mag- 
netic fields, since it is equivalent to the fully polarized 
ferromagnetic state at S = 0, and this polarization gives 
a dominant positive contribution to the Hamiltonian ([2|) . 
Despite being an excited state, the PM phase is stable at 
B ~ [l^ due to spin conservation, which ensures that 
the condensate remains fully polarized. However, after 
introducing magnetic field, the perpendicular part of the 
spin is no longer conserved and the phase becomes un- 
stable, as depicted in Fig.lljb). We note that the growth 
rate of unstable modes of the PM phase at weak fields 
is proportional to the fourth power of the magnetic field 
strength [l9| , and the time of the development of insta- 
bility may be much longer than the condensate lifetime. 

We investigate the destabilization of the PM phase 
in detail by calcula ting the growth rate of linear Bo- 
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+ Sipj with 6ipj = 



goliubov modes 

{uj(t)e^-'' + v*{t)e-'^-'') e~'''^*+"'\ where u{t),v*{t) 
e*"*, w being the eigenfrequency of the excitation. In 
Fig. mc,d) we present spectra of excitations for the 
PM phase in absence and presence of magnetic field. 
The solid line, corresponding to spin quadrupole modes 
[H ESl J exhibits a minimum, which in nonzero magnetic 
field. Fig. [TJd), gives rise to unstable modes (imaginary 
frequencies) localized around a finite wavelength value, 
of the order of the spin healing length . Roton modes 
are characterized by an energy minimum at nonzero mo- 
mentum 0. Here the minimum corresponds to massless, 
gapless (A, p, = 0) rotons at i? = 0. This picture is 
analogous to that at the point of instability in dipolar 
condensates [l3|- However, the phonon part of the spec- 
trum at low momenta and the maxon (maximum) feature 
[l3| are absent in the current case. 

To illustrate the onset of instability in the PM phase, 
in Fig. [2] we present results of numerical simulations of 
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FIG. 2: Onset of the instability and spontaneous periodic 
pattern formation in a ^^Na condensate in a quasi one dimen- 
sional regime, (a) Density of atoms in the m/ = sublevel, 
no{x,t), in a ring geometry with periodic boundary condi- 
tions at a; = ±80/im. (b) Snapshot of the densities in the 
sublevels at t — 75ms. The upper solid line is the total den- 
sity n, and the lower solid and dashed lines correspond to 
ni and no, respectively. The density n_i is very close to ni. 
Other parameters are iV = 3 x lO", uj± ^2ttx 10^ Hz, M = 0, 
B = 0.4G, A^o = 0.17 N. The excitation spectrum of the initial 
state corresponds to Fig. [lid). 
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^■^Na atoms confined in one-dimensional geometry with 
periodic (ring) boundary conditions [l^. The calcula- 
tions were carried out within the truncated Wigner ap- 
proximation with the initial noise filtered to remove high 
spatial frequencies above fcmax — 0.2/im~^. Such a choice 
of initial noise gives a better agreement between theory 
and experiment 2l|. It is clearly visible that the insta- 
bility leads to spontaneous appearance of transient pe- 
riodic patterns before they dissolve into random struc- 
tures. The spatial period of the pattern corresponds to 
the wavelength of the unstable quadrupole modes. We 
note that similar patterns were observed recently in a 
F ~ 2 condensate prepared in an analogous, transversely 
polarized state subject to magnetic field 22]. 

In Fig. [3] we present results in two- and three- 
dimensional geometries with harmonic trapping poten- 
tials. Following a typical experimental scenario, the ini- 
tial state is prepared from the ground state with all the 
atoms in the m/ = 1 component, by coherently trans- 
ferring part of the atoms to the other components and 
keeping the mutual phase equal to 9 — 0. Because the 
density of atoms varies in space and so the atom distri- 
bution in the PM phase, it is not possible to prepare a 
perfect initial state with this technique. We prepare an 
imperfect PM state by setting the atom distribution to 
A''_|_:7Vo:-^-=2:l:2, which corresponds approximately to 
the average density (n). Figures 3(a)-{c) show snapshots 
of the atomic density in the nif = sublevel in the two- 
dimensional setting with increasing condensate dimen- 
sions. The rotonlike instability gives rise to the emer- 
gence of regular patterns in similar conditions as in the 
untrapped case, but only if the condensate size is larger 
than the pattern characteristic length, which roughly cor- 
responds to the healing length for the average density 
(n). The symmetry of the patterns reflects the symme- 
try of unstable modes, which can change substantially 
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FIG. 3: Transient patterns formed spontaneously in a quasi 
two dimensional and three dimensional settings in harmonic 
traps. Pictures (a)-(e) show snapshots of the density and (f) 
shows the isosurface of the density of atoms in the m/ = 
sublevel after a specified time of evolution. The initial states 
are close to the PM state (see text), except (d) and (e), where 
the initial state is the 2C state with wide instability spectrum 
in (d) and narrow, bounded instability spectrum in (e). Pa- 
rameters are (a) = 3 x 10'', ijJ± =2nx 80 Hz, t = 75 ms, (b) 
= 16 X 10*, uj± = 27vx 50 Hz, t = 57ms, (c) TV = 2 x 10^ 
cjx = 27r X 35 Hz, t = 62 ms, (d) as in (c) with t = 48 ms, (e) 
iV = 4 X = 27r X 100 Hz, t = 46ms, (f ) = 3 x 10^ 

LJx = = 2-K X 100 Hz, t — 152 ms. Other parameters are 
B = 0.4G except (e) where B = 0.8G, and oj^ = 2n x 10^ Hz 
in (a)-(e). Densities in trap centers correspond approximately 
to Bo = 0.4 G. The frame (e) is magnified two times. 



depending on the ratio between the characteristic length 
and the condensate size. For comparison, in Fig. ^d) 
we present irregular pattern generated from the 2C state, 
which has a non-roton instability spectrum, with no lower 
bound for wavevectors k. Figure [3Ke) shows a crystalline 
pattern created from the 2C state at a higher magnetic 
field, after the bifurcation with the PM state, when the 
spectrum becomes bounded (see below). In Fig. [31Jf ) we 
present an example of a tetrahedron-like pattern gener- 
ated in a three-dimensional setting. We note that almost 
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FIG. 4: (a) Emergence of the narrow instability spectrum 
(shaded) in the PM phase. At the instability threshold, 
-B = 0, unstable wavevectors are localized around a nonzero 
value fco. (b) For comparison, the instability of the 2C phase 
has initially no lower bound for the wavevector length, and 
produces irregular patterns, see Fig. [Sj^d). At stronger mag- 
netic fields, the instability spectrum becomes bounded and 
gives rise to regular structures as in Fig. |3je). 



low magnetic fields. This argument also allows to under- 
stand the absence of rotonlike instability in ferromagnetic 
condensates. In this case, the excited nematic and 2C 
states [HI are characterized by zero transverse spin, and 
the spin conservation does not prevent energy relaxation. 

In conclusion, we demonstrated that the excitation 
spectrum of an antiferromagnetic condensate in a simple 
model with pure contact interactions can exhibit a ro- 
tonlike minimum. Under the influence of magnetic field, 
this minimum gives rise to an instability and can lead 
to spontaneous emergence of regular periodic, polygonal, 
polyhedral or crystalline patterns. Theoretical considera- 
tions suggest that the appearance of rotonlike instability 
is related to spin conservation, which inhibits the relax- 
ation of excess energy at low magnetic fields. 

This work was supported by the EU project NAME- 
QUAM and the ARC Centre of Excellence ACQAO. 



perfect crystalline patterns could be generated form the 
PM state in the absence of trapping potential. 

The above results can be understood by considering 
the spectra of unstable modes in Fig. 2] Both the in- 
tuitive picture and the presented numerical simulations 
suggest that the regular structures emerge when the Bo- 
goliubov spectrum contains only a limited range of un- 
stable modes, with wavevector lengths close to the char- 
acteristic length of the generated pattern. This kind of 
spectrum appears naturally when the instability is a re- 
sult of crossing of a roton minimum into the unstable 
domain. We argue that this situation occurs in excited 
states whose relaxation is prevented by a conservation 
law. If the constraint can be relaxed gradually, the insta- 
bility is likely to emerge only in a narrow range of finite 
wavevector values. This intuitive reasoning is supported 
by the following calculation of the wavelength of the most 
unstable Bogoliubov mode of the PM state. At the in- 
stability threshold, B = 0, we choose one of the possible 
forms of the mean- field PM phase tpo = ■^), 

and the quadrupole mode (5-0 — £(5; l)^*'^'^ with- 
out loss of generality. To estimate the value of unstable 
wavevectors, we calculate the normalized spin energy per 
atom Cg = ■^C2n, — 8e^C2n + O(e^) and the kinetic en- 
ergy per atom Ck = 2e^fi,^fc^/TO -I- O(e^), assuming that 
a single most unstable linear mode is dominant. Since 
the spin energy is converted into kinetic energy, the en- 
ergy conservation condition h^kg = 4mc2f^ determines 
the value of fco, which agrees perfectly with that of un- 
stable modes appearing at B = in Fig. Ufa). In com- 
parison, the instability of the 2C state, which is initially 
not separated by an energy gap from the ground state, 
cf. Fig.[lja), has no lower bound for the wavevector value, 
see Fig.mjb). Hence the narrow bounded spectrum in the 
PM phase is due to finite energy gap and spin conserva- 
tion, which inhibits relaxation of the excess energy at 
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